Observation of entanglement between two light beams spanning an octave in optical frequency 
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We have experimentally demonstrated how two beams of light separated by an octave in frequency can be- 
come entangled after their interaction in a \^ nonlinear medium. The entangler consisted of a nonlinear 
crystal placed within an optical resonator that was strongly driven by coherent light at the fundamental and 
second-harmonic wavelengths. An inter-conversion between the fields created quantum correlations in the am- 
plitude and phase quadratures, which were measured by two independent homodyne detectors. Analysis of 
the resulting correlation matrix revealed a wavefunction inseparability of 0.74(1) < 1 thereby satisfying the 
criterion of entanglement. 



Central to modern techniques in optical metrology has been 
the ability to make connections between light beams that span 
an octave in optical frequency. This development has real- 
ized the optical-comb whose offset-frequency can be directly 
linked to the SI definition of the second jl]]. As such, spec- 
troscopic measurements can now be made with an absolute 
accuracy beyond one part in 10 15 Q], and have enabled un- 
precedented testing of fundamental quantum mechanical ef- 
fects J^t]. Parallel to this has been progress in the area of quan- 
tum optics, where sources of non-classical light have pushed 
the sensitivity of spatial and temporal interferometric mea- 
surements beyond their associated classical bounds J3l. Re- 
cent achievements in the strength and stability of squeezed 
light sources suggest their imminent application to improving 
the sensitivity of ground-based gravitational wave detectors 
JH. Correspondingly, a source of entangled light with beams 
that span an octave in optical frequency and which share quan- 
tum correlations in amplitude and phase, called harmonic en- 
tanglement, has the potential to be applied to the heterodyne 
stabilization of optical-combs used in metrology @]. 

The goal of producing harmonic entanglement has only re- 
cently been actively pursued by experimental groups world- 
wide. Milestones include demonstrations of two-color 
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FIG. 1 : (color online) Entanglement generation between the reflected 
fundamental and second-harmonic light of an optical parametric am- 
plifier (OPA). 



discrete-variable entanglement from cascaded two-photon 
emission, and spontaneous parametric down-conversion J3l; 
and two-color continuous-variable entanglement spanning 
1 ran from non-degenerate optical parametric oscillators 
(NDOPO) . The large spans offered by harmonic entan- 
glement were predicted to occur in second-order nonlinear 
processes such as second-harmonic generation (SHG), above- 
threshold optical parametric oscillation (OPO), and optical 
parametric amplification (OPA) HElH. Although a quan- 
tum correlation in amplitude between the fundamental and 
second-harmonic fields of SHG has been observed fl2ll . con- 
firmation of entanglement also requires the detection of phase 
quadrature correlations of the (necessarily) bright light beams. 
Techniques have been used to rotate the phase quadrature onto 
the amplitude quadrature using either an under-coupled opti- 
cal resonator, or an unequal arm length Mach-Zehnder inter- 
ferometer 1 1311 . Yet with the most recent results from NDOPO 
appearing limited by excess phase noise 11411 . a demonstration 
of harmonic entanglement has remained elusive. 

In this Letter, we present observations of entanglement in 
continuous-variables between two modes of light at 532 nm 
and 1064 nm spanning an octave. Our best results show a 
wavefunction inseparability of 0.74(1) < 1 between the re- 
flected fundamental and second-harmonic fields of a degener- 
ate OPA. We overcame the problem of excess phase noise by 
choosing the driving field parameters such that the noise effec- 
tively cancelled. In addition, we demonstrated the technique 
of optical carrier rejection (OCR) using impedance-matched 
cavities, which enabled us to directly perform sideband homo- 
dyne detection on light beams that were originally an order 
of magnitude brighter than the local oscillator. From mea- 
surements of the amplitude and phase quadratures we deter- 
mined the correlation matrix elements and tested for entangle- 
ment. We exploited the flexibility of our experimental setup 
to smoothly drive the system from parametric amplification 
to de-amplification. The results are in good agreement with a 
model of pump-depleted OPA. 

The OPA process is described by an efficient phase- 
dependent exchange of energy between a fundamental seed 



2 



OPO J< 



CO. 




whose inverse is found: 



-5 5 

Seed Field Amplitude ai n [VmW] 



FIG. 2: (color online) Map of entanglement (inseparability) as a 
function of seed and pump driving fields as based on our OPA model. 
Symbols mark where regions have been experimentally probed. 



and second-harmonic pump field. Due to this process, the light 
fields reflected from an OPA can be entangled in their ampli- 
tude and phase quadratures. Our aim is to obtain expressions 
for these observables and apply the inseparability criterion to 
verify entanglement. We extend our previous model ll ill to 
include a term that describes the effect of guided acoustic- 
wave Brillouin scattering (GAWBS). Consider a mode of light 
a and its second-harmonic b, (wavelengths A a = 2Ab), which 
interact via a nonlinearity of strength e in a single mode cavity 
with total decay rates K a ,b- The intra-cavity fields are coupled 
to the environment through a mirror K a i,bi, and more weakly 
via other loss mechanisms n a 2.b2- The system is driven by two 
coherent fields a in , f3 in and can be modeled by IU6I1 



a = — (n a +i5w a )a + ea/b + A[ n , 
b = -(n b +i5w b )b - ^ea 2 + B in , 



(1) 

(2) 



where input fields are denoted by A ln = \f2K a \a- m + 
J2 ■y/2K a ,jSA j A n ; B in = y/2n bl (3i n + ^2 ^2Kb,jSBj !iti with 
j€{l, 2}, K a = n a i+K a 2 and K b = Hbi+ ^b2- The GAWBS noise 
is modeled by the detuning terms Sw a j, = (— 2irc/ A ai b)£, a ,bSP 
which are driven by the dimensionless noise term SP hav- 
ing variance one, but coupled via £ a .&; see lfl5tl . We work in 
the Heisenberg picture where the annihilation operators a and 
b (and corresponding creation operators) evolve, from which 
the amplitude and phase quadrature operators are constructed, 
X + =a^ + a and X~=i((v — a), respectively. For simplicity, 
we drop the hat notation henceforth. Operator linearization is 
used to obtain the fluctuations (SX^ , SX^~) centered around 
the classical steady-state solutions (a, (3) Jig]. Fourier trans- 
forming into the frequency domain allows one to solve for the 
driving fields in terms of the intra-cavity fields. This depen- 
dence is reversed when the equations are expressed in a matrix 
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(3) 

where {5X±,5X^} and {SX A in , SXg in } are the intra- 
cavity and accumulated input field quadratures, respec- 
tively; A± = K a — ifl ± e|/3| cos (9/3, B = — e\[3\smdp, 
C = — e\a\ cos6 a , D = — e\a\s'm9 a , E = K b — 
iVt, F a , b = 2i\a,f3\sin6 a ,(3(-2irc/\ a , b )t;a,b, G a , b = 
2\a,f3\ cos Q ^(-2TTc/XaM)£,a,b, with 6 a = Arg(a),^ = 
Arg(/3); and ft the sideband frequency. The fields re- 
flected from the resonator can be directly obtained using the 
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input-output formalism, SX A1 rcf 

The resulting bi-partite Gaussian states are completely 
described by the correlation matrix of elements C^ n = 

\ (SXtSX l n + SX l n SXt), given by 
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where {k,l} £ {+, — } and the the reflected field notation has 
been simplified with {to, n} G {^4l rc f a, -Bl rc f i— > b}. We 
use the quantity of inseparability X as a measure of entangle- 
ment 1 1811 . It is defined from M ab in the following way: 
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The inseparability criterion T<1 is a necessary and sufficient 
condition for entanglement provided that the correlation ma- 
trix has been brought into standard form, which is obtained by 
applying local symplectic transformations r a ,b to each mode 
separately (SX^>— > e ±ra SX a t etc.) to minimize I and satisfy 
k + =k~=k; see jl^l . Inseparability takes values in the range 
[0, oo) with zero signifying perfect entanglement, and unity 
the classical limit. 

A map of inseparability as a function of seed and pump 
fields is plotted in Fig. 2. The entangled states produced by 
the OPA are biased, which means that the correlations of both 
quadratures are not of equal strength. These biased states, 
however, may still be optimally entangled (X— > 0), in the ideal 
limit of a lossless system [20]. Except for the GAWBS coef- 
ficients £ a ,b which were fitted to the phase quadrature mea- 
surements, all parameters used in our model were obtained by 
characterization of the experimental setup 112 111 . 

Figure. 1 shows a schematic of the experiment. Generating 
and measuring harmonic entanglement demands the construc- 
tion of an OPA optimized for the conflicting requirements of 
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FIG. 3: (color online) Inseparability as a function of angle parameter 
with total input power held at 65 mW. Entanglement is achieved 
for values < 1. Solid and dashed curves are from our model with 
and without GAWBS noise, respectively. Measurements marked by 
(circles)crosses are (un)corrected for artifacts from the OCR process. 



high escape efficiency and low threshold-power. Only bright 
driving fields can yield strong entanglement, but the measure- 
ment of phase on bright fields is technically difficult. The 
source for the experiment was a Nd:YAG laser with internal 
frequency doubler, producing continuous-wave light at wave- 
lengths 1064/532 nm. Where "/" refers to parameters for 
each wavelength in that order, and which we call red/green to 
aid readability. The light first needed spatial and temporal fil- 
tering via transmission through an optical cavity of linewidth 
0.4 MHz for red; and two stages of filtering with 1.0 MHz 
then 1.9 MHz for green. Both beams were near shot-noise- 
limited at 7.8 MHz, showing only an excess noise variance of 
0.1% per mW of power. A fraction of the light was diverted 
for the local oscillator beams (LO), before preparing the seed 
and pump beams with a set of RF sidebands for servo control. 

The entangler itself consisted of a temperature controlled 
10 mm long PPKTP crystal placed at the smaller focus of a 
bow-tie cavity geometry designed to minimize astigmatism. 
A dispersion plate placed in the larger focus allowed adjust- 
ment of the red-green phase shift. Characterization of the cav- 
ity yielded finesses of 57/8.8; linewidths 18//60MHz; and 
reflection coefficients 74/70%. The input-output mirror re- 
flectivities were 90/53%; total losses from the other mirrors 
were 0.6/0.9% and losses due to the crystal were 0.2/6.7% 
per round trip; giving an escape efficiency of 92/86%. The 
seed/pump were s-polarized and mode-matched into the cav- 
ity to 98/99%. The cavity length was controlled to resonate 
with the seed field. The pump field was made to co-resonate 
by adjusting the crystal temperature. Turning the dispersion 
plate optimized the effective nonlinearity and lowered the 
OPO threshold pump power to P t h = 85(5) mW. Note that 
our doubly-resonant cavity ensured a low P tn while distribut- 
ing losses fairly equally over both colors. After interaction in 
the cavity, the reflected red and green beams (now entangled) 
were separated, and a 1% tap-off on each was detected. This 
enabled locking of the seed-pump relative phase so that the 
OPA could operate either as an amplifier or de-amplifier. 



The vital step in detection was to remove the carrier light 
while preserving the sidebands. This was achieved by us- 
ing near impedance-matched filter cavities having linewidths 
0.4/1.0 MHz. It was possible to reject 21/26 dB of car- 
rier light, with the cavity length locked onto resonance, 
and the eigenmode of the OPA spatially mode-matched to 
99.9(1)/99.9(1)%. Two homodyne detectors received the 
remaining reflected light which contained the sidebands, but 
also a small contribution of the original carrier light in higher 
order spatial modes originating from mode mismatch of the 
pump and seed into the OPA eigenmode. 

Each red/green homodyne detector had a pair of photodi- 
odes with estimated quantum efficiency 95(2)/95(5)%. The 
LO was mode-matched to the OPA eigenmode with visibility 
99.0(2) /99.4(2). The total detection efficiency was estimated 
to be 87/88% and the shot-noise to dark-noise clearance was 
15 dB. The homodyne LO phase was servo-locked to either 
the phase or amplitude quadratures of the signal beam. Before 
taking entanglement readings, the locked quadrature angles 
were aligned to coincide with the max/min noise power. This 
method guaranteed orthogonality of the quadratures being 
measured, and allowed us to infer that the cross-quadrature 
correlations of a single color (C^~, etc.) were zero. We have 
achieved a homodyne condition with signal to LO power of 
1:>30, representing up to 26 dB filtering of the carrier field 
whilst preserving almost all the sideband signal. 

Two independent electronic channels were prepared, 
through which the photo-currents of the homodyne detectors 
were band-passed, amplified, and mixed down at 7.8 MHz, 
before being sampled at 44.1 kHz to acquire 2.6 x 10 5 points. 
A complete data set contained paired measurements of the ho- 
modyned signal beams (SX^, 8Xj~), the signal beam without 
LO for OCR noise characterization, and also the vacuum noise 
and detector dark noise references. From this set, it was pos- 
sible to normalize the quadrature amplitudes and obtain the 
correlation matrix elements. Corrections were made for the 
offsets incurred by dark noise, and excess-noise from higher- 
order spatial modes in the signal beam (which are artifacts of 
the OCR process, and not the entangler itself). 

Results: Our aim was to drive the entangler across a range 
of parametric processes covering OPO, SHG, and OPA. We 
fixed the total driving field power to 65 mW (which is 76% 
of Pth), and adjusted the balance of power between seed and 
pump. In Fig. 3, wavefunction inseparability is plotted as 
a function of the polar angle of Fig. 2. By moving away 
from pure OPO, entanglement was observed over a broad 
range of angles (— 0.41, +0.15)7T, which covered paramet- 
ric de-amplification through SHG, and approached the neutral 
point (where net inter-conversion of fundamental and second- 
harmonic is zero). The effect of pump-enhancement in this 
region was prominent. Most of the fundamental field was 
converted into second-harmonic field. We also found en- 
tanglement in a narrow range (+0.41, +0.47)7r correspond- 
ing to parametric amplification with weak pump-depletion. 
The entanglement observed in the pump-enhanced region was 
slightly better than in the pump-depleted region with I = 
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FIG. 4: (color online) (a) Time-series quadrature data showing cor- 
relations (scaled 50%). (b) Shaded ellipses follow a contour of the 
resulting probability distribution. Dashed circles mark the quantum 
noise limit. The quantum correlation in amplitude is evident since 
the ellipse falls within the circle, (c) Dual quantum correlations are 
exhibited when the correlation matrix is brought into standard form. 



0.76(2) and I = 0.79(2), respectively. The dashed curve in 
Fig. 3 is a prediction based on a previous model 111 111 , while the 
solid curve is from the extended model with GAWBS noise. 
The latter agrees well with the experimental results. 

We observed the best harmonic entanglement in the para- 
metric de-amplification region of strong pump-enhancement. 
The seed and pump powers were set at 81 mW and 9 mW, re- 
spectively (see * in Fig. 2). We repeated measurements of the 
correlation matrix for eleven runs over many days. The en- 
semble average of those matrix elements in linear scale, and 
their 95% confidence intervals based on the run-to-run vari- 
ability, are presented here: 



M a 



0.71(1) -0.25(1) -0.02(6) 

2.45(12) -0.07(10) +1.42(5) 

-0.25(1) -0.07(10) 0.83(2) 

-0.02(6) +1.42(5) 2.56(6) 



The matrix revealed that both colors were amplitude squeezed 
with C++= 0.71(1) and C++ = 0.83(2). The phase quadra- 
tures showed anti-squeezing of C~~= 2.45(12) and C^~ = 
2.56(6), which imply that the Heisenberg uncertainty rela- 
tion was satisfied well above the minimum uncertainty bound. 
These apparent "mixed" state statistics are a requisite of har- 
monic entanglement. To compute the inseparability, we per- 
formed local symplectic transformations r a = 0.11(1), r& = 
0.15(2) numerically to each mode such that M & was brought 
into the standard form. Applying the inseparability definition 
Eq. (O, returned a value of X = 0.74(1), thus confirming the 
presence of entanglement. 

A visual representation of the correlations within the entan- 
gled state is shown in Fig. 4(a), where time-series quadrature 
data of the second-harmonic field was plotted against the fun- 
damental field. The ellipse in (b) marks the standard deviation 
contour of the resulting probability distribution. The quan- 
tum anti-correlation in amplitude is evident as the ellipse falls 
within the circular boundary set by vacuum states. For the 



phase quadrature, only a classical correlation can be seen, but 
the proximity of the phase correlations to the classical bound 
is sufficient for the preservation of entanglement. This feature 
is symptomatic of biased entanglement flii EoTl . In (c) we 
performed local symplectic transformations to bring the cor- 
relation matrix into standard form. This led to the amplitude 
quadratures becoming correlated and the phase quadratures 
anti-correlated by an equal amount, thereby optimally redis- 
tributing the quantum correlations over both quadratures. 

In conclusion, we have observed entanglement between two 
light beams spanning an octave that was generated from an 
optical parametric amplifier. The results are in good agree- 
ment with a theoretical model that was tested across a range 
of parametric processes. An entanglement source of this type 
could aid the linking together of two optical systems at dif- 
ferent wavelengths. In particular, bridging the octave in op- 
tical frequency between telecommunication windows and the 
prominent multi-level atomic systems used in quantum mem- 
ories and information processing [22]. 
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